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OO ■ We apply dynamical systems techniques to investigate cosmological models in- 
O , 

. spired in scalar-tensor theories written in the Einstein frame. We prove that if the 

O ' potential and the coupling function are sufficiently smooth functions, the scalar field 

Q ' almost always diverges into the past. The dynamics of two important invariant sets 



is investigated in some detail. By assuming some regularity conditions for the poten- 



O ', tial and for the coupling function, it is constructed a dynamical system well suited 

^ ' to investigate the dynamics where the scalar field diverges, i.e. near the initial sin- 



gularity. The critical points therein are investigated and the cosmological solutions 
associated to them are characterized. We find that our system admits scaling solu- 
tions. Some examples are taken from the bibliography to illustrate the major results. 
Also we present asymptotic expansions for the cosmological solutions near the initial 
space-time singularity, which extend in a way previous results of other researchers. 

PACS numbers: 04.20.Jb, 04.20.Dw, 98.80.-k, 98.80.Es, 95.30.Sf, 95.35. +d 



I. INTRODUCTION 

It is believed that the field responsible for early time inflation as well as the field driving 
the current accelerated expansion is in the form of an scalar field (see jl| for the observational 
status of the acceleration of expansion). Theories including scalar fields, such as scalar- 
tensor theories (STT) of gravity [2|, [sj], can be supported by fundamental physical theories 
like superstring theory [4]. Otherwise, scalar fields can be viewed merely as convenient 
heuristic models to elucidates certain (qualitative) dynamical features of the early and/or 
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the late time universe ^. Quintessential dark energy (DE) models 6|, for instance, are 
(heuristical proposals) described by an ordinary scalar field minimally coupled to gravity. A 
wise selection of the scalar field self- interacting potentials can drive the current accelerated 
expansion. Other scalar field models have been treated in the literature (see for instance, 
the reviews j?!). 

The natural generalizations to models of quintessence evolving independently from the 
background matter are models which exhibit non-minimal coupling between both compo- 
nents. Experimental tests in the Solar system impose severe restrictions on the possibility 
of non- minimal coupling between DE and ordinary matter fluids jsj. However, we argue 
that, due to the unknown nature of dark matter (DM), it is possible to have additional 
(non gravitational) interactions between the DE and the DM components. This argument 
does not enter in conflict with the experimental data, but, when the stability of dark energy 
potentials in quintessence models is considered, the dark matter-dark energy coupling might 
be troublesome j^. 

The effective interaction dark energy-dark matter appears when we apply conformal trans- 
formations in scalar-tensor theories (STT) ^. About these theories, it is known that they 

n 

survive several observational tests including Solar System tests [10| and Big-Bang nucleosyn- 
thesis constraints [l^. Its simplest proposal is the Brans-Dicke theory (BDT) Q], in which 
a scalar field, Xi ^cts as the source for the gravitational coupling with a varying Newtonian 
'constant' G ~ X^^- More general STT with a non-constant BD parameter i^(x), and non- 
zero self- interaction potential V^(x), have been formulated, and also survive astrophysical 
tests js]. 

The action for a general class of STT, written in the so-called Einstein frame (EF), is 

n 

given by [131 ]: 



(1) 



In this equation R is the curvature scalar, is the a scalar field, related via conformal 
transformations with the dilaton field, x- ^(0) is the quintessence self-interaction potential, 

^ See 0] and references therein for the analysis of a scalar field responsible for both the early and the late 
time inflationary expansion in the context of brane cosmology. 



^ See the reference 
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for applications of conformal transformations in both relativity and cosmology. 
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is the coupling function, Cmatterif^,^ fJ', x{(p)~^9a/3) is the matter Lagrangian, /x is a 
collective name for the matter degrees of freedom. 

By considering the conformal transformation gaf3 = xi^fy^^Qap and defining the Brans- 
Dicke coupling 'constant' uj{x) iii such way that d(j) = ±^/iu{x) + 'i/'^X^^dx and recalling 
^(x) = the action ([T]) can be written in the Jordan frame (JF) as (see jl^): 



SjF = d^'x^li^^xR - I'^i^xY - V{X) + CmatterifJ^, V/X, (2) 



'Ma 

Both frames are both formally and physically equivalent 



15j. This fact removes previ- 



ous doubts and fully establishes equivalence at the classical level; however, this does not 
guarantee physical equivalence at the quantum level [l6| . 

By making use of the conformal equivalence between the Einstein and Jordan frames we 
can find, for example, that the theory formulated in the EF with the coupling function = 
Xoexp((0 — 0o)/cc), zu = ±^yuJo'+3j2 and potential V{(j)) = (3 exp{{a — 2)zu / {(j) — (po)) 
corresponds to the Brans-Dicke theory (BDT) with a power law potential, i.e., uj{x) = 
^ix) = I^X"'- Exact solutions with exponential couplings and exponential potentials (in 



the EF) were investigated in 17l |. 



In the STT given by ([2]), the energy- momentum of the matter fields is separately con- 
served. However, when written in the EF ([T]), this is no longer the case, although the overall 
energy density is conserved. In fact in the EF we find that 

where 

Ta(3 = --^^^ |v^X"^>C(/i, V/i,x"W)} • 



In the present investigation we study FRW space-times with flat spatial slices, i.e., we 
consider the line element: 

ds^ = -de + a{tf [dr^ + [dO^ + sin' Odip^) ) . (3) 

We use a system of units in which SttG = c = h = 1. 

The simplest way of modelling matter and energy in the Universe is to consider that the 
energy- momentum tensor T^p is in the form of a perfect fluid 

Tp = diag(-p,p,p,p) , 
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where p and p are respectively the isotropic energy density and the isotropic pressure (con- 
sistently with FRW metric, pressure is necessarily isotropic [3]). For simplicity we will 
assume a barotropic equation of state p = (7 — l)p. 

In this case, the vector field = has only one non-null component, the time component, 
so it can be written as Qp = (Qo, 0,0,0). Qo can be considered in some way as a rate of 
energy exchange between the scalar field and the background. 

One of the first papers to take seriously the possibility of interaction in scalar field cos- 
mologies, from the dynamical system perspective, was 19(]. In that paper we can find a 
review on the subject. There it was investigated the interaction terms (in the fiat FRW 
geometry) 6 = —a(f)p and 5 = apH, where a is a constant, (p is the scalar field, p is the 
energy density of background matter and H stands for the Hubble parameter. The first 
choice corresponds to an exponential coupling function x(0) = Xoexp (2a;0/(4 — 87)) . The 
second case corresponds to the choice x = XoO'~'^°'^^'^~^"'^ (and then, p oc a°'~^^), where a 
denotes the scale factor of the Universe (recall that the former derivations are only valid 
in a fiat FRW model). Other phenomenological coupling functions were studied elsewhere. 
We want to draw the attention of the reader to a physically well motivated approach to 



the coupling function in 



20[ |. In that paper it was investigated a coupling term of the form 



S = —ap, where a is a constant (F in their notation). As commented in that reference, if 
a > 0, the model can describe either the decay of dark matter into radiation, the decay of 
the curvaton field into radiation or the decay of dark matter into dark energy (see section 
III of {2^ for more information and for useful references). In the reference 21 1, the authors 
construct a family of viable scalar-tensor models of dark energy (which includes pure f{R) 
theories and quintessence). They consider a coupling between the scalar field and the non- 
relativistic matter in the Einstein frame of the type -in our notation- x(0) = e"^*^*^, with 
Q constant. By investigating a phase space the authors obtain that the model posses a 
phase of late-time acceleration preceded by a standard matter era, while at the same time 
satisfying the local gravity constraints (LGC). In fact, by studying the evolution of matter 
density perturbations and employing them, the authors place bounds on the coupling of the 
order \Q\ < 2.5 x 10~^ (for the massless case). By a chameleon mechanism the authors show 
that these models can be made compatible with LGC even when \Q\ is of the order of unity 
if the scalar-field potential is chosen to have a sufficiently large mass in the high-curvature 
regions. 



About the dynamics of coupled dark energy, it was found (see jl^ and references therein) 
that typically at early times (t —>■ 0) the BDT solutions are approximated by the vacuum 
solutions and at late times {t —>■ oo) by matter dominated solutions, in which the matter is 
dominated by the BD scalar field (denoted by x the Jordan frame). Exact perfect fluid 
solutions in STT of gravity with a non-constant BD parameter uj{x) have been obtained 



by various authors (see [221). Coupled quintessence was investigated also in [23|] by using 
dynamical systems techniques. 

In order to classify the global behavior of the solutions of it is required a detailed 
knowledge of the form of the scalar field potential (and of the coupling function x) ■ However, 
up to the present, there exist no consensus about the specific functional form of V{(f)) (and 
of x(0))- As a consequence it would be of interest to classify the dynamical behavior of 
solutions without the prior of specifying the functional form of the potential function (and 
of the coupling function). In the literature on General Relativity (GR) several attempts 
have been made in this more general direction. The dynamical systems techniques have 



proven very useful to do so 



24 



25 



26!]. In this respect we know the dynamical behavior of 
scalar field space-times for a wide class of non-negative potentials (within Einstein gravity 
(EG)) \ 

In this investigation we want to study, from the dynamical systems point of view, a 
phenomenological model inspired in a STT with action ([1]), where the matter and the 
(quintessence) scalar field are coupled in the action ([1]) through the scalar tensor metric 



x{4')~^9ai3 131]. ^ We consider arbitrary functional forms for the self-interaction potential 
and the coupling function for the scalar field (p. When we take the conformal transformation 
allowing writing the action in the JF as in ([2]) the coupling function x should be interpreted 
as the dilation (BD) field and the corresponding uj{x) as the varying BD parameter. 

By employing Hubble-normalized dynamical variables in addition to the scalar field we 
find, by investigating the phase space, that the scalar field almost always diverges into the 
past, allowing us to identify this regime with the physical region on a vicinity of the initial 
space-time singularity. For finite values of the scalar field we find that the late time (early 



In the reference 



27|, the nonnegativity of the potential is relaxed and new results within the context of 



EG generalizing those in [2q are obtained. 
^ In [A] it is provided a summary of the main results of the theory of dynamical systems that are used in 
this paper. 
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time) attractor is associated with the minimum of the logarithm of the potential (coupling) 
function. 

By assuming some general regularity conditions for the potential and for the coupling 
function when oo, and using the formalism developed in 2^, we are able to explore 
the phase space corresponding to this limit. We are able to compute the critical points 
corresponding to that limit and to characterize the cosmological solutions associated to 
them. Scaling solutions do arise in this regime. Also we are able to characterize the initial 
singularity in STT. In|A]an example taken from literature is given to illustrate the formalism 
developed for analyzing the region where the scalar fields diverge. 

II. THE MODEL 

In the model described by the action ([T]), the background energy density does not nec- 
essarily corresponds to a Dark Matter component and, analogously, the scalar field energy 
density does not necessarily corresponds to a Dark Energy component. However, as done 
many times in the literature (see for instance, 17|, l28|), we assume that is a quintessence 
scalar field, which is coupled metrically to a background of a perfect fluid. This model is 
viable phenomenologically. The possibility of a universal coupling of dark energy to all sorts 



of matter, including baryons (but excluding radiation) is studied in [29l |. 

A. The field equations 

By using the line element ([3]), the Einstein's field equations (derived by varying the action 
([I])) are: a) the Raychaudhuri equation 

H = -l{lP + ^'), (4) 

where the dot denotes derivative with respect the cosmic time t, p is the energy density of 
dark matter, b) the Friedmann equation 



3H' = ^^' + V{<P)+p, 

and the continuity equation 



(5) 
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p = -37//p-^ (4-37)p0^. (6) 
This equation can be integrated in quadratures to give the useful relation 

The equation of motion of the scalar field (written as two differential equations, one for 
and the other for 6) reads 



where the coma denotes the derivative with respect to the scalar field. 

Observe that, if x = const., the equations for the minimally coupled theory are recovered. 

From the equations (jH El ED we see that {H, p, 0, 0) G remains in the hypersurface 
defined by the restriction ([5]). Thus, defining an autonomous system in the phase space 

n=\{H, p, 0, 0) e : = ^02 + v{<P) + p| (9) 

We shall consider the following general assumptions: V{(f)) G and V{(j)) > 0, ^ C*^ 
with x(0) > 0, p > 0, < 7 < 2, 7 7^ 4/3 instead of considering specific choices for both 
the potential and the coupling function. 



III. QUALITATIVE ANALYSIS ON THE HUBBLE NORMALIZED STATE 

SPACE 

In this section we rewrite equations (jH O [TJ [8]) as an autonomous system defined on a 
state space by introducing Hubble-normalized variables. This variables satisfy an inequality 
arising from the Friedmann equation (Q. We analyse the cosmological model by investigating 
the flow of the autonomous system in a phase space by using dynamical systems tools. In 
order to make the paper self-contained we offer in |X] some terminology and results from the 
theory of dynamical systems that we use in the demonstration of our main results. 
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A. Normalized variables 

In order to analyze the initial singularity and the late time behavior it is convenient to 
normalize the variables, since in the vicinity of a hypothetical initial singularity physical 



variables would typically diverge, whereas at late times they commonly vanish 30 1. 
Let us introduce the following normalized variables 

and the time coordinate 

dr = 3Hdt. (11) 

Besides, if no additional information is available on the functional forms of the coupling and 
the potential, the most natural variable to add to the former ones is the scalar field itself. 

B. The autonomous system 

The field equations (jH El El E]) can be used to obtain evolution equations for the variables 
lOl) and the scalar field (h. 



x' = h (2y2 + z^^) , (12) 

s 1^2 N x'^dAyU) (4-37)^2 
y' = y' + l^ [z'l -2)y- + ^^^^a^lnx(0), (13) 

z' = \z {2y- + {z^ - 1) 7) - ^i^^9,lnx(0), (14) 

^y (15) 

where the prime denotes derivative with respect to r. This is an autonomous system where 
the variables are subject to the constraint 

y^ + z'^ + l/3xV(0) = 1. (16) 

Observe that y"^ + z"^ < 1, since V{(j)) is nonnegative and the restriction (fT6ll holds. 
By the |2] in |A] it is possible to prove that any combination of the sets 

x<0, x = 0, x>0,z<0, z = 0, z>0 
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is an invariant set of the flow of our dynamical system provided x is at least of class C^. 
Taking into account this result we can restrict our attention to the flow restricted to the 
phase space: 



x,y, z) eM} : x>Q, z>Q,y'^ + + l/3a;V(0) = l} . (17) 

Observe that if x 7^ and V{(f)) > we can use the constraint f|T6l) as a deflnition for x. 
In this way the evolution equation of x decouples from the evolution equation of the other 
variables. The equation ( fT3l) can be rewritten as 



'^\lny(0) + ii-^9,lnx(0)+y^ + ^(A-2)y. (18) 



Hence, we may then study the dynamical system on given by the equations (fT2!) . (fT8 
and (|T4|) with the constraint + z^ < 1. 



C. Main Lemma 

First we want to prove the lemma [Tl 1 . 1 1 st at ing that the orbits passing through an arbitrary 
point p G interpolates between a regime where the Hubble parameter diverges 

(containing an initial singularity into the past) to a regime where the background density is 
negligible (into the future). Those orbits represent cosmological solutions with non- vanishing 
dimensionless background energy density and flnite and positive Hubble parameter. The 
result stated by the above lemma is obtained by constructing a monotonic function deflned 
on an invariant set and by applying the LaSalle monotonicity principle (Theorem [2] in Rl) . 
This lemma is very important as a tool for investigating the past attractor or cu-limit set of 
the flow (see deflnitions [T] and [2] in El) : it is necessarily located at the invariant set {x = 0}. 

Lema III.l Let be = {{x,y, z,(j)) E : z > 0} and let be X'^ = {{x,y, z,(j)) E H : x = 
0}. Then for all p G — X^ the a- and uj-limit sets of p are such that a{p) C X^ and 
uj{p) C where dZ^ denotes the boundary of Z^ . 

Proof. By the proposition [2] in 1X1 we have that S = Z^ — X° = {(x, y, 0) G S : 2; > 
0,a; > 0} is an invariant set of the flow (I12til5p . Let the function 

^(x,l/,^,0)=(^)'x(0)-^+^-/^ 
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with > 0, be defined on S. It is a monotonic decreasing function (in tlie direction of the 
flow) in S, since its directional derivative through the flow is Z' = — 7Z (which is obviously 
negative). ^ The rank of Z is (0, 00). Let he s ^ S — S = dZ^ U X^. It is verified that 
Z{s) — s> as s — s> dZ^ and Z{s) ^ 00 as s — > X^. Hence, applying the LaSalle monotonicity 
principle (see Theorem O in the E}, we have for all p G Z~^ — X° that a{p) C and 
uj{p) C dZ~^ , as required. □ 



D. The dynamics restricted to the invariant sets x = and z = 0. 

Now we will to characterize the dynamics in both the invariant sets = {{x, y, z, 0) G 
S : X = 0} and = {(x, z, 0) G S : z = 0}. 

The dynamics in the invariant set X^ is governed by the differential equations 

.^i(.-.^)(..-.,.^i^^), 



3?/, (20) 



plus the algebraic equation: 



y' + z^ = 1. (21) 

The only critical point of this system (with (p bounded) is the critical point Q with 
coordinates (x, y, z, (p) = (0, 0, 1, (pi) with x'{4>i) = and x{4>i) 0- 

The critical point Q represents matter dominated cosmological solutions with the Hubble 
parameter diverging. Since d^xi^Pi) = 0? they are solutions with minimally coupled scalar 
field (with negligible kinetic energy). The potential function is also unimportant in the 
dynamics. 

The eigenvalues of the liberalization around Q are |, 7, Ai ± y^A^ + A2 '^(^^^^ ; where 
Ai = (-2 + 7)/4 < 0, and A2 = (4 - 87) /6. 

Then, the local stability of Q in the invariant set X° is as follows (we are assuming that 
the barotropic index 7 satisfies < 7 < 2): 

1. Q is a stable focus if < 7 < 4/3 and x"(0i) < -^^^ or 4/3 < 7 < 2 and 
X"(0i) > 



^ See definition 2] in the \M 
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2. g is a stable node if < 7 < 4/3 and > x"{<Pi) > or 4/3 < 7 < 2 and 
< X"ih) < 

3. Q is a saddle point if < 7 < 4/3 and x"{4>i) > or 4/3 < 7 < 2 and x"{4>i) < 0. 

4. Q is non hyperbolic, if = 0, in which case, there exists a 1-dimensional stable 
manifold which is tangent to the axis y at Q (the associated eigenvector is By = (1, 0)). 
There exists also a 1-dimensional center manifold tangent to the line (1 — 7/2)?/ — 

v/2730 = at g. 

The dynamics in the invariant set is governed by the differential equations: 



and ( l20l) . plus the equation 

+ l/3xV(0) = 1, (23) 

where V{(f)) is given as input. 

The only critical point (with bounded) in the invariant set Z° is the critical point 
P with coordinates x = J v{<f>2) ' V — ^1 ^ — ^1 — 4>2 with x{4>2) 7^ 0, V^'(02) = 0. The 



eigenvalues of the linearization around P are: 0, — |, — | ±1^1 — |V^"(02)/^(02)- The zero 
eigenvalue has the associated eigendirection e^,. 

The local behavior of the critical point P in the invariant set is as follows: 

1. P is a saddle if V"{(j)2) < 0, 

2. P is a stable node if < V"{(t)2) < |V'(02), and 

3. P is a stable focus if V"{(p2) > |^(02)- 

4. P is non hyperbolic (in the invariant set Z^) ii V"{(f)2) = 0- 

When the orbits located at Z^, approach this critical point, the energy density of DM and 
the kinetic energy density of DE tends to zero. In this case the energy density of the Universe 
will be dominated by the potential energy of DE. Hence, the Universe would be expanding 
forever in a de Sitter phase. 
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E. Conditions for the divergence of the scalar field backwards in time 

Now we will prove the Theorem IIII.2I which, essentially, states that if the potential and 
the coupling function are sufficiently smooth functions, then for almost all the points lying 
in a 4-dimensional state space, the scalar field diverges when the orbit through p is followed 
backwards in time. This theorem is an extension of the theorem 1 in the reference 2J] to 
STT. 

Theorem III. 2 Assume that ^^^^^ ^(0) '^'"^ positive functions of class C^. Let p be a 
point in S, and let 0~{p) he the past orbit of p under the flow of ( I7M73]) with constraint 
^W\) . Then, is almost always unbounded on 0^{p) for almost all p. 

Proof. Let p G S such that (f) is bounded on 0~{p). It follows that 0~{p) is contained on 
a compact subset of (the closure of) S. Hence, that trajectory must asymptotically approach 
some limit set a{p) (see the analogous of theorem [T] in El for a limit sets). From f|T2|) follows 
that a; is a monotonic increasing function through the fiow and then, it must be constant in 
a{p). 

There are two possibilities i) ?/ = 2; = at a(p) or ii) a; = at a{p). 

The only invariant set with y = z = (and bounded) is the critical point P 

Now we will prove that P is not the past asymptote of an open set of orbits of S. 

Observe that at least one of its associated eigenvalues has always negative real part. 
Hence, by the Center Manifold Theorem (theorem [1] in [A]), we can conclude that there 
exists an invariant stable manifold, E"^, of P, intersecting P. The existence of an stable 
manifold of dimension r > implies that all solutions asymptotically approaching P in the 
past (i.e., those ones approaching P as r ^ —00) are in an invariant unstable manifold or 
a center manifold of dimension 4 — r < 4. 

Then, the only possibility is that a{p) is contained in x = (i.e., it is contained in the 
circumference y"^ + z"^ = 1.) 

The only invariant sets for the fiow given by fll9ti2ip are the hypersurfaces |?/| = 1 and 
the critical point Q. Then, there are two possibilities: a{p) = Q or the a-limit set of p lies 
on the hypersurface \y\ = 1. In the last case, from ( l20l) . we have that is unbounded on 
a{p), a contradiction. 

In order to complete the proof we need only to demonstrate that Q is not the past 
asymptote to an open set of trajectories in S. 
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In fact, at least one of its associated eigenvalues has negative real part. In view of the 
Center Manifold Theorem, we can conclude that there exists a local stable manifold, E^, 
that intersects Q of dimension s > 0. All orbits close to Q in this set exponentially approach 
Q as r +00. As before, the existence of a local stable manifold of dimension s > (equal 
1 or 2) implies that all solutions past asymptotic to Q must lie on an unstable manifold or 
center manifold of dimension 4 — s < 4. □ 

Theorem IIII.2I allow us to conclude that in order to investigate the generic asymptotic 
behavior of the system fll2m5p with restriction (fl6D it is necessary to study the region where 
= ±00. However, as has been investigated in j32] (where results from 25I are extended), the 
region = ±00 is not exclusively associated to the asymptotic behavior towards the past. In 
fact, the scalar field can diverge towards the future, provided additional requirements under 
the potential and under the coupling function are fullfilled. But, this will be the purpose of 
a forthcoming paper. 

F. The flow near cj) = +00 

In this section we will investigate the flow near = +00 following the nomenclature and 



formalism introduced in 24] ■ Analogous results hold near = — 00. 

By assuming that V,x ^ with exponential orders and M (the set of all class k WBI 
functions; see the definitions [5] and [6] in [B]) respectively, we can define a dynamical system 
well suited to investigate the dynamics near the initial singularity. We will investigate the 
critical points therein. Particularly those representing scaling solutions and one associated 
with the initial singularity. 

Let C S be the set of points in S for which > e^^, where e is any positive constant 
which is chosen sufficiently small so as to avoid any points where or x = 0, thereby 
ensuring that Wv{'^) and W^{ip) are well-defined (see definitions of W's and of hatted 
functions in lBj) . 

We now make the coordinate transformation 



(x, y, z, 0) '^-^^ (x, y, z, ^) (24) 



on S^, where /(0) tends to zero as tends to +00 and has been chosen so that the conditions 
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i)-iii) of definition [6] in [B] are satisfied with k = 2. 

Substituting this new coordinates in the equations f|T2|) . f|T8|) and f|T^ we obtain the 
3- dimensional dynamical system: 

y' = y' + l{z'l-2)y- ^'-y^''\ Wy + N) + '-^^^^{W, + M), (25) 

= Iz {2y^ + {z^ - 1) 7) + '-^^^^ [W, + M) , (26) 

^' = \j\7y- (27) 

We may identify with its projection into so that we have = 

{0 < if < f{e~^), < y"^ + z"^ < 1} . The variable x can be treated as a function on defined 
by the constraint equation which becomes 



y'' + z'' + l/3x^V{ip) = 1. (28) 

The directional derivative of x along the flow generated by (12511271) may be obtained directly 
by equation (fT2|) . 

Since /', Wy and are at = we may extend (12511271) onto the boundary of to 
obtain a system on the closure of S^, i.e., S^. From definition [6], /', Wy and W^^ vanish 
at the origin and are each of second order or higher in (p and /' is negative on E^. 



1. Critical points. 

The system (I25H27I) admits the critical points labelled by Pi, i G {1,2,3,4,5,6}. In the 
following we discuss the existence and the stability conditions for the critical points. In the 
table [I are displayed the values of some cosmological magnitudes of interest for the critical 
points (the deceleration parameter, the effective equation of stare (EoS) parameter for the 
total matter, etc). 

1. The critical point Pi with coordinates y = —1, z = and (f = exists for all the 
values of the free parameters. The eigenvalues of the linearized system around Pi are 

Ai_i = 2 — a/2/3A^, Ai^2 = — ^"^^2^5^^^ ^^'^ ~ ^' ^^^^^ critical point is 
non hyperbolic, then the Hartman-Grobman theorem does not apply. By the Center 
Manifold Theorem there exists: 
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(a) an stable invariant subspace of dimension two (tangent to the y-z plane) if: i) the 
potential is a WBI function (see definition [5] in [B]) of exponential order > 
and the coupling function is a WBI function of exponential order M < — ^^T^^^ 
(provided < 7 < |), or ii) the barotropic index satisfies | < 7 < 2, the potential 
is a WBI function of exponential order N > \/6 and the coupling function is a 
WBI function of exponential order M > — ^^^-"4^'' i 

(b) an unstable invariant subspace of dimension two (tangent to the y-z plane) pro- 
vided the potential is a WBI function of exponential order < -\/6 and the cou- 
pling function is a WBI function of exponential order M such that M > — ^^'^"4^'* 
(respectively, M < — ^['^'^^ ) provided < 7 < | (respectively, | < 7 < 2); 

(c) a 1-dimensional center manifold which is tangent to the critical point in the 
direction of the axis ip. This center manifold can be 2-dimensional or even 3- 
dimensional (see the discussion on the point 3). 

2. The critical point P2 with coordinates y = 1, z = and if = exists for all the 
values of the free parameters. The eigenvalues of the linearized system around P2 are 



^2,1 = 2 + a/2/3A^, and A2,2 = Ai,2 and A2,3 = (see point 1). Hence the critical point 
is non hyperbolic, then the Hartman-Grobman theorem does not apply. By the Center 
Manifold Theorem there exists: 



(a) an stable invariant subspace of dimension two (tangent to the y-z plane) if: i) 

/6(7-2 
37-4 



< -VQ, M > for < 7 < |, or ii) I < 7 < 2, < -VQ and 



(b) an unstable invariant subspace of dimension two (tangent to the y-z plane) pro- 
vided A^ > -^6, and M such that M < (respectively M > 
provided < 7 < | (respectively | < 7 < 2); 

(c) a 1-dimensional center manifold which is tangent to the critical point in the 
direction of the axis (p. This center manifold can be 2-dimensional or even 3- 
dimensional (see the discussion on the point 3). 

In the following section we shall study the initial spacetime (big bang) singularity. 
The critical points Pi^2 can account for that singularity. They are in the same phase 
portrait for the values -^6 < A^ < ^6 and -^^|^ < M < ^^f^ and < 7 < 
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(in which case they have a 2-dimensional unstable manifold and a 1-dimensional center 
respectively). It is easy to show that the Hubble parameter (and the matter density) 
of the cosmological solutions associated to these points diverges into the past. The 
scalar field also diverges, it equals to +00 (respectively —00) for Pi (respectively P2). 
However, even in this case, the past attractor corresponds to Pi since /' < and for 
y > the orbits enter the phase portrait and P2 acts as a saddle. The last point can 
be a past attractor only on a set of measure zero (if (p — 0). 



3. The critical point P3 with coordinates y = '^(^2+7) ' ^ ~ ^(^2+7)^ (p — 

exists if < 7 < I and -4^^ < M < 4^^. The eigenvalues of the ma- 
trix of derivatives evaluated at the critical point are A3^i = ^^^~^''i2(-^2)^~^^^ ' ^^3.2 = 
+ (M + N)M + + 7, and Ag.g = 0. Hence the critical point is non 

hyperbolic, then the Hartman-Grobman theorem does not apply. Under the above 
existence conditions we find, by the Center Manifold Theorem, that there exists a 
stable manifold of dimension two for the values of the parameters: i) M < and 
N > or ii) M > and N < ^^2^37-1)""'^' " Otherwise there 

exists an unstable manifold of dimension one (in this case the stable subspace is 1- 
dimensional) . The center manifold is in both cases 1-dimensional. If M = T^^q^^ 
this critical point reduces to Pi 2- In this case the center subspace is 2-dimensional 


1 



and is spanned by the eigenvectors 



VV 







. The center manifold 



is tangent to the center subspace at the critical point. If additionally |A^| = -s/6, the 
center manifold is 3-dimensional. 

4. The critical point P4 with coordinates y = z = and (p = exists if |A^| < \/6. 

Observe that this point reduces to Pi_2 if = 6. The matrix of derivatives evaluated 
at the critical point has the eigenvalues A4,i = | {N'^ - 6) < 0, A4,2 = lN{2M + N) - 
j{MN + 2)^ and A4^3 = 0. Hence the critical point is non hyperbolic and, as before, the 
Hartman-Grobman theorem does not apply. However, we can use the Center Manifold 
Theorem to investigate the stability of this critical point. The structure of the center 
manifold is as follows: 

(a) if A4^i < and A4_2 7^ the center manifold is spanned by e^. Then, it is 1- 
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dimensional. Before analyzing this case in detail, we will provide additional 
information about the structure of the center manifold. 

(b) if M = ^^^3^_^4)'^ and N"^ < 6, the center subspace is spanned by the eigenvectors 



ez and e,^. 



(c) if A^^ = 6 and M ^ T-^^z^q^^, it is spanned by the eigenvectors e.y — 



-4+37 

and e,^. 





VV 



(d) if A^^ = 6 and M — ^-^^z^q^^, the center manifold is 3-dimensional. 

The local behavior described in the cases above (excluding the first case) is in some way 
special. It requires fine tuning of the free parameter. However, the typical behavior 
(in the invariant set 2; = 0) is the existence of a one dimensional center manifold 
Cjv through P4, which is tangent to the z-axis (if A4^i < and ^ 0). Cjv is 
an exponential attractor on a sufficiently small neighborhood of P4. It is intuitively 
obvious from the geometry (for instance, observe figure 1), that any solutions past 
asymptotic to P4 must lie on the center manifold. 

Let us investigate the case in which A4^i < and A4^2 7^ 0. Of course, in this case 
the stable manifold is at least 1-dimensional (and as we mentioned before the center 
manifold is 1-dimensional). 

The structure of the stable subspace is as follows: 

(a) if the potential is of exponential order zero [N — Q), then, the critical point has 
coordinates (0,0,0). The eigenvalues of the linearization are (—1,0, — |) and in 
this case, the stable subspace is generated by the eigenvectors 62^,62; 



(b) if < 7 < I, -76 < < 0, and M > 



2(Af2-37) 
Ar(37-4) ' 

2(Ar2_3^) 



or 



(c) if I < 7 < 2, -v^ < < 0, and M < -j^j^; or 

(d) if < 7 < |, < A^ < |, and M < 4/^^; or 

(e) if I < 7 < 2, < A?" < Ve, and M > '^l^^^.^J-l the stable subspace is generated 
by the eigenvectors 6^,62. 
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(f) By interchanging > and < in the inequahties for M in the cases (b)-(e) we find 
that the stable manifold is 1-dimensional and is tangent to the critical point in 
the direction of e^^ (accordingly, the unstable subspace is spanned by e^). 



5. The critical points Ps^e with coordinates y = m{3'^4)-2N ^ ^ = =F '^^^%^+M\'i-3'!)^ 
(respectively) exists if the following conditions are simultaneously satisfied: 4N{2M + 

N) - 6{MN + 2)7 >0,T {2N + M(4 - 87)) > and ^^'^(SMci-s'T))^''^' ^ ^ (i"^-' 
the critical points are real- valued, and they are inside the cylinder S^). 

The associated eigenvalues are 



A± - " + 
^5,6 - — ± 



^8 (/?2 + 2772) a2 _ 2/5(7 - 4) (/52 - 2I672) a - (7 - 2) {f3^ - 2X672)^ 



f3 6^6/57 

and A5,6 = 0, where a = 3 (A^(7 - 2) + A/(37 - 4)) and /3 = 2 (2A^ - M(37 - 4)) . 
Assuming that the conditions for existence are satisfied, we can analyze the stability 
of the critical points by means of the Center Manifold Theorem. We find that the 
non null eigenvalues can not be either complex conjugates with positive real parts or 
real- valued with different sign, then the unstable subspace of P^^q is the empty set. 
Then, the stable subspace is 2-dimensional (provided A5 6 is the only null eigenvalue). 
When the orbits are restricted to this invariant set, the point P^^q acts as an stable 
spiral (if the eigenvalues are complex conjugated) or as a node (if the eigenvalues are 
negative reals). The conditions on the parameters for those cases are very complicated 
to display here. 

If M = "^^^^^^J-l the Ps^e reduces to P4 but in this case, the center manifold is 2- 
dimensional and is spanned by e^, e^,. 



G. The flow near (j) = —00 

With the purpose of complementing the global analysis of the system it is necessary 
investigate its behavior near = —00. It is an easy task since the system (jU El [3, [H]) is 
invariant under the transformation of coordinates 
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Point 


y 


Z 






Acceleration? 


Pi 


-1 





1 


1 


no 


P2 


1 





1 


1 


no 


P3 


5 




6^ 


7 + (7 -1)5 


2 

< 7 < — and 


|M| < r 


Pa 


N 

VE 





1 


^ 3 




P5,6 


6\/67 




2{2a+l3) , 4327 , 1 


(7+2)/32+4a(7+l)/3-43272 


/3 ^ 3 


13 




/37 + -32- + i 


/3^7 



TABLE I: The properties of the critical points for the system ()25ll27p . We use the notations 

M(37 - 4) 

^(7 - 2) 



a = 3 {Ni^ - 2) + M(37 - 4)) , /3 = 2 (2iV - M(37 - 4)) , <5 = ^^^^ ^^^^ , and F = . 



(0,0) ^-(0,0), y^f/, x^S, 

where t/(0) = V{—(j)) and S(0) = x(— 0). Hence, for a particular potential V, and a partic- 
ular coupling function x, the behavior of the solutions of the equations (jH O [71 [8]) around 
= — 00 is equivalent (except for the sign of 0) to the behavior of the system near (p = 00 
with potential and coupling functions U and S, respectively. 

If U and H are of class S'^, the preceding analysis in can be applied (with and adequate 
choice of e). 

In the following we will denote to the set of class functions well behaved in both 
+00 and —00. We will use Latin uppercase letters with subscripts +00 and —00, respectively 
to indicate the exponential order of functions in +00 and in —00. 
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H. The global geometric structure of the phase space 

Let be fl{xo) the region of the phase space given by (fT71) with x < Xq, then, since x is 
monotonic decreasing, this set equals the union of its past orbits. 

The procedure to define a coordinate system near — oo, given in the section IIII Gl can 
be used to embed fi(xo) as a compact differentiable 4-dimensional manifold S(xo) such that 
the vector field defined by fll2til5p can be smoothly extended over S(xo). 

With this purpose we define an atlas as follows. 

First, the interior of S(xo) is defined as the set 

{{(j),x,y, z) G : < x < Xo,y^ + z"^ < 1 } 

where we use the local chart (coordinate system) given by ( fTOl) . Obviously, this set is 
bounded in the variables x, y, z. We define a second local chart {ip, x, y, z) in the open 
subset of this set for which x > e~^, for e small enough, by fl2^ . An analogous local chart 
{ip, X, y, z) can be defined near = — oo by the procedure given in the section IIII Gl 

The construction is completed by attaching a boundary which is defined taking the union 
of a; = 0, a; = Xo, = and the circumference + 2;^ = 1 to each local chart. 

By construction, S(xo) is compact and it is embedded in M"^. 

The vector field defined by ( fT2lfT5l) can be smoothly extended over the boundary of 
S(xo) such that S(xo) is the union of its past orbits, i^(xo) is a 3-dimensional hypersurface 
embedded in S(xo). It is important to note that ^{xq) approaches the non-physical boundary 
along the intersection of the plane x = with the plane = and the circumference 
y"^ + z"^ = 1. This set is called non-physical boundary of and denoted by dVt. 

IV. THE INITIAL SPACE-TIME SINGULARITY 

In this section we will study the initial space-time (Big-Bang) singularity. The critical 
points Pi 2 can represent such a singularity. They live at the same phase space for the values 
of M, N and 7 in the intervals -^6 < < ^6, - ^^1"/^ < M < ^^1"/^ and < 7 < | 
(in this case, they have an unstable 2-dimensional manifold and a center 1-dimensional 
manifold). It is easy to show that the Hubble parameter and the matter energy density of 
the associated cosmological solutions diverge towards the past. The scalar field diverges too, 
and it is equal to +00 and to —00 for Pi and P2 respectively. However, even in this case. 
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the possible past attractor corresponds to Pi since /' < whereas for ?/ > the orbits enter 
the phase space and P2 acts as a saddle. The critical point P2 can act as a past attractor 
only in a set of initial conditions of measure zero (when = 0). 



A. Analysis near Pi. 

From the analysis in section [Tll F 11 it seems reasonable to think that the initial space-time 
singularity can be associated to the critical point Pi. Its unstable manifold is 2-dimensional 
provided < Vg. The asymptotic behavior of neighboring solutions to Pi can be approxi- 
mated, for r negative large enough, as 

y(r) = -1 + 0(e^^'i"), z{t) = 0(e^i-^"). (29) 

By substitution of (129|) in (fTSll . and integrating the resulting equation, we obtain 

Hr) = ^{-r + 4>)+0{e'^-n- (30) 
Then, by expanding around r = — cxd up to first order, we get 
2 / ~\ ^.1.9^ A . f [2 

" ' ' ^ I I r^f „M 



where h denotes higher order terms to be discarded. 

Then we have a first order solution to (I25ti27p . Also, by substitution of ( l29ll in (fT2l) and 
solving the resulting differential equation with initial condition x(0) = xq we get the first 
order solution 

X = xqc^ . (31) 

Then, we have t — ti = ^ f x{T)dT = l/3xoe'^. For simplicity let us set ti = 0. 
Neglecting the error terms, we have the following expressions 



(32) 



where c = I/Sxqc'^. This asymptotic solution corresponds to the exact solution of dUElElIH]) 
when V vanishes identically and x is a constant (the minimal coupling case). Hence, there 
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exists a generic class of massless minimally coupled scalar field cosmologies in a vicinity of 
the initial space-time singularity. 

The above idea can be stated, more precisely, as the 

Theorem IV. 1 Let be V G S'^ with exponential order N satisfying N < -\/6 and let be 
X ^ £+ with exponential order M such that 

i) 0<-f<l andM> -^^^^ or 



l<^<2 and M< -^^^^ 



Then, there exist a neighborhood JV (Pi) of Pi such that for each p G A/'(Pi) the orbit ipp past 
asymptotic to Pi and the associated cosmological solution is: 

H = 1 + O (evit)) , (33) 



3t 

lhni + 0{tev{t)), (34) 
^t-' + Oievit)), (35) 



'--2 
2 ^ 



p = xlt-'x{-Jlln^A il + Oitevit))) (36) 



where ey(t) =tV (^-^Ini 



Before proceeding to the proof of this theorem, let us make a few comments. 



Since V & has exponential order A^, then, by applying Theorem 2 in 



24| we have 



limrV^(-W-ln-) = lim e-^^^VU) = 0, Va > \ -N. 
t^o V 3 c (p^oo V 3 

Then, for N < \^ the error terms O {ev{t)) and O (tey(t)) are dominated by the first 



order terms. If < ^ / | both error terms tend uniformly to zero. 



2 

~ 37 
_ In i 1 

3 



On the other hand, since x ^ has exponential order M, then, t ^x(—-i/| In- 



tends uniformly to zero, as t — > in the cases i) < 7 < I and M > or ii) | < 7 < 2 

and M < In the above cases, the matter tends uniformly to zero as Pi is approached. 
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a. Proof of Theorem IIV.II From the equation ([T2l) and using (l28i) as a definition 
for y we get the equation: 

Using the first order expressions z = 0{e^^'^'^) and x = x^e^ we have the differential equation 
for x: 



d\nX 1 277/ N 2r , , 

Where we denote by h any collection of higher order terms to be discarded. 
By integrating in both sides of fl37j) we find the solution 



(37) 



X = xqc exp 
XqC I 1 



Xn 



V{if)e^^ 1 + h 



3Ai,i 



^0 t7^,„N„2t 



3Ai,i 



Viip)e'^ + h 



In the above deduction we have used the auxiliary result proved in |24| : 



(38) 



All 



(39) 



which is valid if G with exponential order N < We have used, also, the approxi- 
mation e" ~ 1 + M. 

Now we want to derive an second order expression for t. 



^F(^)e^-) + h. 



3 " V 9Ai 



In the above deduction we have used the auxiliary result proved in 2J]: 



V{^)e^^dT = V(<^)e^" + h 



The equation poj) may be inverted, to second order, to give 



(40) 



(41) 



xqc'^ = 3 i t + 



Ai,i 



Substituting this result in (l38i) we get: 



and then. 



x{t) =3t + h. 

^ ' x{t) 3t 3Ai,i 



The equation fl25|) can be written as 

dlny / 7\ 2 1 2tt/ x (1 - Z/^ - -2^^) ^ttt ,a 2:2(4-37).— 



By the same arguments as in the deductions of (!38|) and (H2|) we get 



™2 
„2t 



-1 + 7TT^e^-r(y.) + h 
-^1,1 



Combining the expansions (jH]) and (H2|) in 0(t) = we find 



0(t): 

This equation can be integrated up to second order to get 



The equation (l26l) can be written as 



dhiz 1 n— , , ,^ 7,,^ n, ■»(— 4 + 37)— , , 
— = --xV(^) + (1 - I)(l - + l^L_^,y,(^,, 

where we have used the constraint equation (128|) as a definition for y"^. 
By the same arguments as in the deductions of (|38|) and P2|) we get 



By definition 

w / N x'(/~Hv?)) , , 
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Then, by using the first order expression 



r\v) = <P= \j\{-r + 0o) + O(e^M-) + h 

(as derived in former sections) and integrating out the resulting expression with respect to 
r we get the estimation 

j WMdr = -^\nx{v)dT -Mr + h. (48) 



By substitution of (HHll in ( 1471) we get 

Z = ZqC 



Xot'-hi^)"^-' + h (49) 



where xo = zo {f-) ^^"^ • 

Combining the expansions for z and x in p = ^ we find 



P = ^ (l + ^^('^)^') x(0)^"'xoV^ + (50) 
Observe that the second term, h, on the right-hand of fH6|) tends to zero when t 0. 

I In 4 

3 c 



This allows to Taylor expand V and x around 0* = —^/| In 4 to get 



V{(j){t)) = V^(0*)(1 + aWv{(l)*)V{(f))f) + h (51) 

and 

xm)) = xma + aw^{cpnv{<p)t') + h (52) 

where a is a constant. By substituting equations fl^ and ([S2D into the equations fll5l 113, 
H6l [50]) the theorem is proven. □ 

B. A global singularity theorem 
Finally, we will state (without a rigorous proof) a global singularity theorem which is in 



some way an extension of Theorem 6 in 2J] (page 3501). It is not totally an extension of 
this theorem, since in our framework it is very difficult to prove that the correspondence 
with the massless minimally coupled scalar field cosmologies is one-to-one. 
The theorem states the following: 
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Theorem IV. 2 Let be V e S'^ such that < 6 and x ^ such that 



< 7 < I and M±^ > 



V6(7-2) 
37-4 



or 



a) I < 7 < 2 and < 



^6(7-2) 
37-4 



Then, it is verified asymptotically that: 



where ey{t) = tV (± 




(53) 



(54) 



(55) 



(56) 



Sketch of the proof 




Following the reasoning [24'], it is sufficient to demonstrate that almost all the solutions 
are past asymptotic to the critical point Pi (in oo or in — oo). Since x is monotonic, it 
is sufficient to consider solutions in Q{xo) C S(a;o) where Xq is arbitrary. Since S(xo) is 
compact and it contains its past orbits, then all the points p must have an a-limit set, 
a{p). Particularly, for the points in the physical space Q{xo), the theorem IIII.2I implies that 
a{p) must contain almost always a critical point with = (0 = icxo). By the discussion 
in the section IIII Ht each point with ip = being a limit point of the physical trajectory, 
must be part of the non-physical boundary dQ{xo) and then must have x = 0. Since x is 
monotonically increasing, the set a{p) must be contained completely in the plane x = 0, or 
namely in dQ{xo). It can be proved that the only conceivable generic past attractor are the 
critical points Pi in ±oo (the other critical points cannot be generic sources by our previous 
linear analysis). 



In this paper we have investigated models with additional (non-gravitational) interaction 
between DE and DM. This kind of interaction is justified if the interacting components are 



V. CONCLUSIONS 
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of unknown nature, as it is the case for the DM and the DE, the dominant components in the 
cosmic fluid. We have investigated these models from the dynamical systems viewpoint. The 
functional form of the potential and the coupling function is arbitrary from the beginning. 
Some general results are obtained and proved by considering general hypotheses on these 
input functions. 

We have proved (by using normalized variables) that the scalar field typically diverges 
into the past. This is formulated in Theorem III1.2I It is an extension of the theorem 1 in 
the reference 2J] to the non-minimally coupled scalar field setting. 



In lemma IIII.ll it is proved that the orbit passing through an arbitrary point p G S 
(representing cosmological solutions with non-vanishing dimensionless background energy 
density and positive finite Hubble parameter) is past asymptotic to a regime where the 
Hubble parameter diverges containing an initial singularity into the past, and is future 
asymptotic to a regime where the background density is negligible into the future. This 
result is obtained by constructing a monotonic function defined on an invariant set and by 
applying the LaSalle monotonocity principle (theorem 4.12, 3l|). 

We have proved the Theorem IIII.2I that makes clear that in order to investigate the 
generic past asymptotic behavior of our system we must seek on the limit where the scalar 
field diverges. 

By assuming some regularity conditions on the potential and on the coupling function in 
that regime we have constructed a dynamical system (well suited to investigate the dynamics 
where the scalar field diverges, i.e. near the initial singularity). The critical points therein 
are investigated and the cosmological solutions associated to them are characterized. We 
find the existence of three critical points P3, P5 and Pq. They are in the boundary of the 
phase space S^. They represent cosmological scaling solutions (where the contribution of 
the dimensionless potential energy is negligible). By tuning the free parameters they can 
be accelerating. In contrast in the reference 2J] there exists only one (in our notation, P4) 
representing an accelerating cosmology. The solutions associated to Pi, 2 (p=f notation 



m 



24| ) represent stiff and then decelerating solutions (actually solutions associated to a 



massless scalar field). 

We have proved a theorem (theorem IIV.1|) which is an extension of the theorem 4 in 24 1 
to the STT framework. Also, we sketch the proof of the global singularity theorem IIV.2[ 
Theorem IIV.2I indicates that the past asymptotic structure of non-minimally coupled scalar 
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field theories with FRW metric, as in the FRW general relativistic case, is independent of the 
exact details of the potential and/or the details of the background matter and the coupling 
function. 

This is a conjecture with solid theoretical and numerical foundations (see figures 1 and 2 in 
IB 21) . To prove that the family of solutions which asymptotically approach Pi are completely 
characterized by the solution space of the massless scalar field cosmological model (i.e., V 
and X s-iid then p, being dynamically insignificant in the neighborhood of the singularity 
Pi) it is required to prove that this correspondence is one-to-one and continuous, which is 
hard to do in our scenario. 
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APPENDIX A: SOME TERMINOLOGY AND RESULTS FROM THE 
DYNAMICAL SYSTEMS THEORY 

In the following x) denotes the flow generated by the vector field (or differential 
equation) 
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x{t) = f{x{t)), x{t) e W. (Al) 

1. Limit sets 

Definition 1 (definition 8.1.1, [36] page 104) A point xq E M" is called an w-limit 
point of X G M", denoted uj{x), if there exists a sequence {ti}, — > oo such that (f){ti,x) —>■ xq. 
a-limits are defined similarly by taking a sequence {ti}, ti — oo. 

Definition 2 (definition 8.1.2, [36] page 105) The set of all w-limit points of a flow or 
map is called a w-limit set. The a-limit is similarly defined. 



Proposition 1 (proposition 8.1.3 j36f |. page 105) Let (pti') be a flow generated by a 
vector field and let M be a positively invariant compact set for this flow (see definition 3.0.3 
page 28 [36]). Then forp G M, we have 

1. i) uj{p) ^ 
a) uj{p) is closed 

Hi) uj{p) is invariant under the flow, i.e., uj{p) is a union of orbits, 
iv) uj{p) is connected. 

A similar result follows for a-limit sets provided the hypothesis of the proposition are 
satisfied for the time reversed flow (see proposition 1.1.14 in [37]). 

2. Center Manifolds 

A general vector field can be transformed locally in the neighborhood of a fixed point 
into a vector field of the form 



x = Ax + f{x,y), (A2) 
y = By + gix, y), (x, y) G x M^ (A3) 



where 
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/(O, 0) = (7(0, 0) = Df{0, 0) = Dg{0, 0) = 0, (A4) 

A is a cx c matrix having eigenvalues with zero real parts, B is a s x s matrix with negative 
real parts, and / and g are functions (r > 2). 

The Center Manifold for the vector field (lASp is defined as follows: 



Definition 3 (definition 18.1.1, 36|] page 246) An invariant manifold will be called a 



center manifold for M3j) if it can locally be represented as follows 

WioM {(x, y) G X R'\y = h{x), \x\ < S, h{0) = 0, Dh{0) = 0} 
for 6 sufficiently small. 

The conditions h{0) = 0, Dh{0) = imply that W{'^^{0) is tangent to E"" at {x, y) = (0, 0), 
where E'^ is the center subspace, i.e., the invariant set spanned by the eigenvectors whose 
associated eigenvalues have zero real parts. 



Theorem 1 (Center Manifold (theorem 2.7.1 38[|)) Let be (pt the flow of a vector 
field, X, then, there exists locally a center manifold, W^^^, containing the origin and in- 
variant under (pt such that Wf^^, has tangent space E'^ at x = 0. This manifold is for all 
A; G N, but its domain of definition can depend on k. Furthermore, there are locally smooth 
stable and unstable manifolds, W^^^ and W^^^, which contain x = 0, are invariant under 0j, 
have tangent spaces E^ and E"^, respectively, and are such that (ptlw^i^ is a contraction while 
(ptlwi^^^ is an expansion. 

3. Monotone functions and Monotonicity Principle 



Definition 4 (definition 4.8 |3ll |. page 93) Let (pt be a flow on M", let S be an invariant 
set of (pt o-nd let Z : S ^ be a continuous function. Z is monotonic decreasing (increasing) 
function for the flow (pt means that for all x & S, Z{(pt{x)) is a monotonic decreasing 
(increasing) function oft. 



Proposition 2 (Proposition 4.1, 3l|], page 92) Consider a differential equation x' 



f{x), X with flow (pt. Let Z : ^ M 6e a (M") function which satisfies Z' = aZ, 
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where a : M" M is a continuous function. Then, the subsets of M" defined by Z > 0, 
Z = 0, or Z < are invariant sets for (f)f 



Theorem 2 (Monotonicity Principle, [311], page 103) Let (pt be a flow on with S 
an invariant set. Let Z : S ^ W be a (M") function whose range is the interval (a, b) 
where a G M U {— oo}, 6 G M U {+cx3}, and a < b. If Z is decreasing on orbits in S, then for 
all X & S, u!{x) C {s G S — S\limy^sZ{y) ^ b} and a{x) C {s G — S\limy^sZ{y) ^ a}. 

APPENDIX B: REGULARITY CONDITIONS OF THE POTENTIAL AND OF 
THE COUPLING FUNCTION AT INFINITY 

In this section we present rigorous statements of what we mean by regularity conditions. 
We present some worked examples that we shall use to construct a toy model. Our purpose 
is illustrated by the techniques for analysis in the region = +oo. In principle the analysis 
is general enough to be applied to more physically interesting situations. 

Definition 5 (see reference [24]) Let : R ^ M 5e a non-negative function. Let 
there exist some 00 > for which V{(j)) > for all (p > (pQ and some number N such that 
the function Wy '■ [(po, oo) R, 

satisfies 

lim Wv{(j)) = 0. (Bl) 
Then we say that V is Well Behaved at Infinity (WBI) of exponential order N . 

It is important to point out that may be 0, or even negative. Indeed the class of 
WBI functions of order is of particular interest, containing all non-negative polynomials 



as remarked in 



24|. 



In that reference it was defined a procedure for classifying the smoothness of WBI func- 
tions at infinity. If we have some coordinate transformation ip = f{(j)) which maps a neigh- 
borhood of infinity to a neighborhood of the origin, then if is a function of (p, 'g is the 
function of whose domain is the range of / plus the origin, which takes the values; 
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lim<^^oofi'(0) , = 



Definition 6 (see reference [24]) A function V is class k WBI if it is WBI and if 

there exists 0o > and a coordinate transformation if = f{(j)) which maps the interval [0o, oo) 
onto (0, e], where e = /(0o) one? lim<^^oo / = 0, with the following additional properties: 

i) f is C^"*"^ and strictly decreasing, 
a) the functions Wv{<^) and f'i^p) are on the 
closed interval [0,e]. 

n 

We designate the set of all class k WBI functions In the table 1 of [24] is incorporated 
a diverse range of qualitative behavior within the framework of £i_ potentials. 



1. Worked examples. 

a. Power-law coupling function. 



Let us consider the coupling function 



= ( ^ ) Xo(</'-</'o)% a > O,const., 00 > 0. (B2) 



Observe that 



dx a{(f)-(j)o) 

for all finite value of 0. Then, the critical point Q located at the invariant set X° does 
not exist for this choice of coupling function. By our analysis in section llllt the early time 
dynamics is associated to the limit where the scalar field diverges. 
This choice produces a coupling BD parameter given by 

Mx) + 3 = ^« f ^ 
3 VXo 



This types of power law couplings were investigated in 331] from the astrophysical view 



point. For scalar-tensor theories without potential, the cosmo logical solutions for the matter 
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domination era (in a Robertson- Walker metric) are a{t) oc {\nt)^" i)/3a^3^ (Int)^. 
The values of the parameter a in concordance with the predictions of are a = 1, 0.33, 3 



(see table 4.2 in 



33|). 



b. The Albrecht-Skordis potential. 



Albrecht and Skordis 3J] have proposed a particularly attractive model of quintessence. 



It is driven by a potential which introduces a small minimum to the exponential potential: 

y(0) = e-^'^(A + (0-5)2). (B3) 

Unlike previous quintessence models, late-time acceleration is achieved without fine tuning 
of the initial conditions. The authors argue that such potentials arise naturally in the low- 
energy limit of M-theory. The constant parameters, A and B, in the potential take values of 
order 1 in Planck units, so there is also no fine tuning of the potential (we suppose also that 
/X 7^ 0). They show that, regardless of the initial conditions, scales, with p oc oc 
during the radiation and matter eras, but leads to permanent vacuum domination and 
accelerated expansion after a time which can be close to the present. 

The extremes of of the potential (IB3P are located at 0^ = ^^^^ ^^"^^ • They are real if 
1 > fi'^A. The local minimum (respectively, local maximum) is located at (respectively 
0"^) since 

±V"(0±) = -2VoVl - Ap2e-(^+^^^v/^^) < 0. 

As we investigated in section IIIIl the late time dynamics in the invariant set is asso- 
ciated with the extremes of the potential. When we restrict ourselves to this invariant set, 
we find that the critical point associated to 0+ is always a saddle point of the corresponding 
phase portrait. The critical point associated to 0^ could be either a stable node or a stable 
spiral if 

(3 + 4^2)2 - ^2 

or 

8(3 + 2p2) 



(3 + V)2 ■ 

The early time dynamics in this invariant set corresponds to the limit = -|-oo. 

Let us concentrate now in how to apply the mathematics introduced in this section. 
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Observe first that 

W^{4>) = = -T-^-^-^ ^ ^xi<P) = (B4) 

a[(p ~ (pQ) (p^+oo 

and 

= d^Vi<l>) /¥{<!>) + /i = ^^'f ~ ^ lim = 0. (B5) 

In other words, the couphng function flB2p and the potential (IBSp are WBI of exponential 
orders M = and = — /i, respectively. 

It is easy to prove that Power-law coupling and the Albrecht-Skordis potential are at least 
under the admissible coordinate transformation 

(^ = 0-i = /(0). (B6) 

Using the above coordinate transformation we find 



2ip 



> 



and 



, ^ = 

2y(Bi/p-l) 

Wviif) = { '^V^^I^^ ' ^ > ^ (B8) 

, (^ = 



fi^) = (B9) 
, if = 



2. The dynamics in the limit (j) +oo. An example. 

We will present here a toy model by choosing a Power-law coupling and Albrecht-Skordis 
potential with the only purpose of illustrate our previous results. Due to the generality of 
our study it can be applied for more realistic scalar-tensor cosmologies. 

Let us consider a toy model with the coupling function (lB2p and the potential (IBSp . 
Observe that this toy model corresponds to a theory written in the Jordan frame with 
action 

SjF= [ d'^Xy^l llxR-l'^i^xf - V{X) + 'C^aUer{f^,Vf^,ga/3)] (BIO) 
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FIG. 1: Orbits in the invariant set {z = 0} C for the model with couphng function ()B2p potential 
()B3p . We select the values of the parameters: e = 1.00, /i = 2.00, A = 0.50, a = 0.33, S = 0.5, and 
(j)Q = 0. Observe that i) almost all the orbits are past asymptotic to Pi; ii) P2 is a saddle, and iii) 
the center manifold of P4 attracts all the orbits in the {z = 0}. However, it is no more an attractor 
in the invariant set 2; > 0, (/9 = (see figure 2 in section lB2]l . 



with BD coupling 



and potential 



where 



( \ 3 2/^ 

^ X =-0 + 0" — 

2 3 VXo 

nx) = xV'^*W(A + (0(x)-i?f) 




3a Vxo 



We have considered also, conformal metric 



gap = ^gap- 
In this example, the evolution equations for y, z, and are given by the equations 
with M = 0, N = -fi and W^, Wy, and f, given respectively by ([BT]), dHH]) and dHl]). The 
state space is defined by = {{y, z, ip) : < y"^ + z"^ < 1, < ip < e} . 

The critical points of the system fl25tl27p in this example are Pi 2 = (t15 0,0), P3 = 
(0,1,0), and P4 = (^,0,0), and P5 



'5,6 = 1 \/|^^ =F— ^2^!"^^^ • points Pi,2,3 exist 
for all the values of the free parameters. The critical point P4 exists for /i^ < 6. The 



critical point P5 exists if yU < — v37 whereas the critical point Pe exists if /i > v37. We will 
characterize the critical points Ps^g in more detail (for the analysis of the other critical points 
we submit the reader to table [T]). The critical points P5 6 corresponds to those studied in the 
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FIG. 2: Orbits in the invariant set {i^ = 0} C for the model with couphng function (|B2p potential 
()B3p . We select the values of the parameters: e = 1.00, /i = 2.00, A = 0.50, a = 0.33, B = 0.5, and 
00 = 0. In the figure i) Pi^2 are local past attractors, but Pi is the global past attractor; ii) ^3^4 
are saddles, and iii) P5 is a local future attractor. 



book (see equation 4.23 p 49) with the identifications ^ = y and = |^ = ^^f^^"'^ 
and k = —fi. As stated in that reference the scalar field 'inherits' the equation of state of the 
fluid, i.e., = 7. Then this solutions represent cosmological kinetic-matter scaling solutions 
^ (the potencial energy density is negligible). These critical points represent accelerating 
cosmologies for < 7 < |. The eigenvalues of the matrix of derivatives evaluated at P^^q are 
^0, — ± ^a/(2 — 7) (247^ + yu2(2 — 97)) j . The orbits initially in the stable subspace of 
Ps^e spiral-in around Ps^e if /^^ > 247^/(— 2 + 97) provided | < 7 < 2,7 7^ |. Otherwise P^^, 
looks like a stable node for the orbits lying in the stable subspace. The center subspace is 
tangent to the critical points in the direction of the ip axis. 



® See reference 



35[ for a notion of 'scaling' solutions, particularly, kinetic-matter scaling solutions. 



